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Abstract 

We study the time evolution of an atom suddenly coupled to a thermal radiation field. As a simplified model 
of the atom-electromagnetic field system we use a system composed by a harmonic oscillator linearly coupled to 
a scalar field in the framework of the recently introduced dressed coordinates and dressed states. We show that 
the time evolution of the thermal expectation values for the occupation number operators depend exclusively on 
the probabilities associated with the emission and absorption of field quanta. In particular, the time evolution 
of the number operator associated with the atom is given in terms of the probability of remaining in the first 
excited state and the decay probabilities from this state by emission of field quanta of frequencies ut- Also, it 
is showed that independent of the initial state of the atom, it thermalizes with the thermal radiation field in a 
time scale of the order of the inverse coupling constant. 

PACS Number(s): 03.65Yz, 05.70.Ln, 05.30. Jp 



1 Introduction 

The study of systems out of thermal equilibrium has been since long time ago one of the main active areas in 
physics. The actual interest ranging from condensed mater physics to cosmology. In most cases the interest is in 
the thermalization process, the determination of the relevant time scales involved, together with an understanding 
of the generation of entropy and particle production in non equilibrium dissipative systems interacting with an 
environment. However, despite the importance associate to these processes, non equilibrium problems are still 
poorly understood [TJ. The nontrivial non equilibrium dynamics of fields, for instance, have diverse applications, 
finding use e.g. in the studies concerning the recent experiments in ultra-relativistic heavy-ion collision [2 ; appli- 
cations to the current problems of parametric resonance and particle production in cosmology [3] ; or in the context 
of the recent studies involving the intrinsic dissipative nature of interacting fields [H [6]. In addition to that, 
typical problems we have in mind to study are those related to the nontrivial out-of-thermal equilibrium dynamics 
associated with phase transitions in different physical systems. As a few examples we may cite include the current 
applications to the study of formation of Bose-Einstcin condensates after a temperature quench [7] , or in the study 
of the dynamics of coupled fields displaced from their ground states as determined by their free energy densities 
[8]. For recent attempts to solve some related problems to the study of systems out of thermal equilibrium see 
Refs. [HI QH IH1 HH EH 031 HE! El HE] , where use has been made of either analytical or numerical approaches in the 
context of specific or general models. For example, numerical studies have been performed in specific field theoret- 
ical models in Refs. [§1 EH [HI Q~2] , where the problems of equilibration and thermalization have been studied. On 
the other hand in Ref. [T3] the role of chaos as a mechanism for quantum thermalization has been considered. By 
supposing the validity of Berry's conjecture [15) it has been showed that a rarified hard-spheres gas approaches a 
Maxwell-Boltzmann, Bose-Einstein or Fermi-Dirac distribution according on wether the wave functions are taken 
to be non-symmetric, completely symmetric or completely antisymmetric functions of the particle position. 

In recent works, in analogy with the renormalized fields in quantum field theory, the concepts of dressed 
coordinates and dressed states have been introduced [TH1 [201 121] • These concepts have been introduced in the 
context of an atom, approximated by an harmonic oscillator, linearly coupled to a scalar field, the whole system 
being confined in a spherical cavity of diameter L. In terms of dressed coordinates, dressed states have been defined 
as the physically measurable states. The dressed states having the physical correct property of stability of the 
oscillator (atom) ground state in the absence of field quanta (the quantum vacuum) . For a recent clear explanation 
see Ref. [25J- Also, the formalism showed to have the technical advantage of allowing an exact computation of the 
probabilities associated with the different oscillator (atom) radiation processes [55]. For example, we obtained easily 
the probability of the atom to decay spontaneously from the first excited state to the ground state for arbitrary 
coupling constant, weak or strong and for arbitrary cavity size. For weak coupling constant and in the continuum 
limit L — * oo we obtained the old know result: e~ rt |19) . Also, considering a cavity of sufficiently small radius 
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)20j . the method accounted for, the experimentally observed, inhibition of the spontaneous decaying processes of 
the atom [221 In Refs. [2J] [23] the concept of dressed coordinates and states have been extended to the case 
in which nonlinear interactions between the oscillator and the field modes are taken into account. Furthermore, in 
Ref. [27] we considered the oscillator electromagnetic field interaction model and in Ref. [28] dressed coordinates 
and states have been introduced in the path integral formalism. 

The aim of the present work is to study the thermalization process in the framework of the afore mentioned 
dressed coordinates and states. The physical situation that we have in mind is an atom (approximated by an 
harmonic oscillator) initially in an arbitrary state, suddenly coupled to a thermal radiation field (approximated by 
an infinite set of harmonic oscillators at thermal equilibrium). Then, the purpose is to study the time evolution 
of this initial state. Fundamental questions that we have to solve are: the atom reaches a final equilibrium state? 
and if this is the case, what is the meaning of this final equilibrium state? In relation with these questions, it is 
also important to know the time necessary for the atom to thermalize with the thermal radiation field. By solving 
exactly this model we expect to gain some insight to solve more complicated problems. 

As already stated, in this paper we will treat a specific model for an atom coupled to a thermal radiation field. 
Initially the system is described by a density operator of the form 



where po is the density operator for the atom, that can be in an arbitrary pure or mixed state and pp is the density 
operator for the radiation field at thermal equilibrium at some given temperature We specify below the form 
of pp. At some time, that we take as t = 0, the atom is suddenly coupled to the thermal radiation field, afterwards 
(the density operator of ) the total system evolves according to the Liouville-Von Neumann equation. An equivalent 
description is to maintain constant the density operator and take the operators (related to the physical observables) 
as time dependent. Then, these operators evolve in time according to the Heisenberg equation of motion 



atom-electromagnetic field system. As a model for this system we consider the one with Hamiltonian given by 



where the limit N — > oo is understood, the subscript refers to the atom approximated by an harmonic oscillator of 
frequency u>o, and k = 1, 2, ...N refer to the harmonic field modes. Also we take uik — 27r/L, Ck = r/uik, V = \[2gKui 
and Aw = u>k+i — LOk = 27r/L, where g is a frequency dimensional coupling constant. At the end we will take the 
continuum limit L — > oo. The last term in Eq. ([3]) assures the positiveness of the Hamiltonian and it can be seen 
as a frequency renormalization of the harmonic oscillator [29l [30] . 

A similar model to the one given by Eq. ([3]) has been used repeatedly from time to time as a simplified model 
to describe the quantum Brownian motion [3TJ [351 1331 133] > the decoherence problem and other related problems 
[351 135] , However, in all these previous works no use has been made of the dressed coordinates. As we explain in 
next section, when considering the Hamiltonian given by Eq. ([3]) as the one for an atom-field electromagnetic field 
system, the introduction of dressed (renormalized) coordinates will be necessary in order to guarantee the stability 
of the atom ground state in the absence of field quanta. 

Along this paper we use natural units h — c = ks = 1. 

2 Dressed (renormalized) coordinates and the dressed density opera- 
tor 

To make this paper self contained in this section we define what has been called dressed coordinates and dressed 
states in Refs. [THl HOI [21] . To understand the necessity of introducing dressed coordinates in the system atom- 
electromagnetic field system described by Hamiltonian ([3]), take = 0. In this case the resulting free Hamiltonian 
admits the following eigenfunctions, 



p = Po ® Pp 



(1) 






(?) 



(g|n ,ni, ...,n N ) 





The physical meaning of ipn n 1 ...n N (q) in this case is clear, it represents the atom in its no-th excited level and rik 
photons of frequencies u>k- Now, consider the state 4'n o...o(q)' the excited atom in the quantum vacuum. We know 
from experience that any excited level of the atom is unstable. The explanation of this fact is that the atom is not 
isolated from interacting with the quantum electromagnetic field. This interaction in our model is given by the 
linear coupling of qo with q^. Obviously, when we take into account this interaction any state of the type V , n o...o(?) 
is rendered unstable. But there is a problem, the state - 0oo...o(g , ) ) that represents the atom in its ground state and 
no photons, is also unstable contradicting the experimental fact of the stability of the atom ground state. What is 
wrong? The first thing that cames in our mind is to think that the model given by Eq. (J3J) is wrong. Certainly, 
we know that the correct theory to describe this physical system is quantum electrodynamics. On the other hand 
such a description could be extremely complicated. If we aim to maintain the model as simple as possible and 
still insist in describing it by the Hamiltonian given in Eq. @ what we can do in order to take into account the 
stability of the atom ground state? The answer lies in the spirit of the renormalization program in quantum field 
theory: the coordinates that appear in the Hamiltonian are not the physical ones, they are bare coordinates. 
We introduce dressed (or renormalized) coordinates, q' and q' k , respectively for the dressed atom and the dressed 
photons. We define these coordinates as the physically meaningful ones. In terms of these coordinates we define 
the dressed states by 



(q'\n ,ni, ...,n N ) d 
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1/4 



(5) 



where the subscript d means dressed state. The dressed states given by Eq. ([5]) are defined as the physically 
measurable states and describe in general, the physical atom in the no-th excited level and rih physical photons 
of frequencies tOk- Obviously, in the limit in which the coupling constant Ck vanishes the renormalized coordinates 
q'^ must approach the bare coordinates q^. Now, in order to relate the bare and dressed coordinates we have to 
use the physical requirement of stability of the dressed ground state. The dressed ground state will be stable if it 
is defined as eigenfunction of the interacting Hamiltonian given by Eq. ([3]). Also the dressed ground state must 
be the one of minimum energy, that is, it must be defined as being identical (or proportional) to the ground state 
eigenfunction of the interacting Hamiltonian. From this definition, one can construct the dressed coordinates in 
terms of the bare ones. Then, the first step in order to obtain the dressed coordinates is to solve for the ground state 
eigenfunction of the Hamiltonian given in Eq. This bilinear Hamiltonian can be diagonalized by introducing 
normal coordinates and momenta Q r and P r , 



N N 
r=0 r=0 

where {tV\ is an orthonormal matrix whose elements are given by [37], 
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i + v ^ 
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with f2 r being the normal frequencies corresponding to the collective modes of the coupled system and given as 
solutions of the equation 



JY 



-fc"r 



In terms of normal coordinates and momenta the Hamiltonian given by Eq. ([3]) reads as 

N 



(8) 



(9) 



then, the eigenfunctions of the Hamiltonian are given by 
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where the subscript c means collective state. Now, using the definition of the dressed coordinates: Voo...o(?') °c 
0oo... o(Q) an d using Eqs. §5§ and (fit))) we get e~ 3 5Zm=o _ e -iE™ from which the dressed coordinates 

are obtained as 



JY 



j - V /— r 



(ii) 



2.1 The dressed density operator 

If no use is made of the dressed coordinates and states, the density operator for the radiation field at thermal 
equilibrium in Eq. ([1]) would be given by 



p fJ = Z 1 exp 



JV 



k=l 



-7. [ a\a k + - 



where dk and a k are annihilation and creation operators and given by 



(12) 



/2uj 



(13) 
(14) 



In Eq. (fT2l) Zp = Ylk=i z 8> * s ^ ne partition function of the thermal radiation field, where 



Zp = Tr fc 



2 sinh 



(^) 



(15) 



Also, the density operator po for the atom would be written in terms of the coordinates qo ■ 

However, as explained above, in the context of an atom-electromagnetic field system and described by Hamil- 
tonian given by Eq. ([3]) it is necessary to redefine what the physical coordinates are for the atom and field modes. 
Then, instead of the density operator given by Eq. (fT2"|) , we have to consider the one written in terms of dressed 
coordinates q' k , as the physically density operator for the radiation field at thermal equilibrium, 



fa = Z R 1 ex P 



-(i^uJk ( a'la' k + 



k=l 



(16) 



where a' k and a' k ^ are dressed annihilation and creation operators and given in terms of the dressed coordinates q' k 
by 




where in position representation p'^ — ~i-^r- Also, the density operator for the atom must be taken as the one 
written in terms of the dressed coordinate q' . 

Now, we are ready to study the time evolution of thermal expectation values for relevant physical operators. 
We will be mainly interested in the present work in the study of the time evolution of the thermal expectation value 
of the time dependent number occupation operator associated with the dressed oscillator (the atom) a' '(t)a' (t). 



3 The thermalization process 

We state the thermalization problem as follows:?,) the initial state given by Eq. (|17|) will evolve in time to a final 
equilibrium state? and ii) if the system evolves to a final equilibrium state, is this an state of thermal equilibrium?. 
Also we would like to know the mean time necessary for the system to reach a final thermal equilibrium state. 

Since any operator can be written in terms of annihilation and creation operators, it will be sufficient to solve for 
the time dependent annihilation and creation operators in order to solve the out of thermal equilibrium problem. 
Using the Heisenberg equation of motion, Eq. we have for the time dependent annihilation operator a'{t), 



^(t)=i[H,a'^t)} (19) 

and a similar equation for a'^(t). Obviously at t — 0, a^(0), is given by Eq. (117j) . This equation can be written, 
using Eqs. ([6|) and (fTTj) . as 

a'(o) = J2 iJ^f^) . (20) 



, jm r L v 2 



In order to solve Eq. (fTi?]) we write a'^(t) as 



iV 



u=0 



«mW = E( b w^ + ^w^ , (2i) 



where -B(i)^„ is a time dependent c-number and the dot means derivative with respect to time. Replacing Eqs. 
([3]) and Eq. (|2Tj) in Eq. (jT9|l . working the commutators and identifying identical operators in both sides of the 
resultant equation, we obtain the following coupled equations for B^ v (t) 



N 9 \ N 

ci 



B^ (t) + |^ + ^-i j B^,(t) - £c fc fl Mfc (i) = (22) 

and 



fc=i fc / fc=i 



S Mfe (t) + w^ft) - CkB^t) = . (23) 

Note that above equations are identical to the classical equations of motion for the bare coordinates that can 
be obtained using the Hamilton equations of motion for the Hamiltonian given by Eq. ([3]). Then we can decouple 
Eqs. (|2"2"|) and (|2"3")) with the same matrix {t r } that diagonalizes the Hamiltonian ([3]), that is, we can write for 

N 

B^^E^fi) (24) 
r=0 

and replacing the above equation in Eqs. (f2"2")l and (|2"3"|) , these equations decouple into 

cj(t) + n?cj;(t) = o ) (25) 

from which we obtain C£(i) = a^e^* + fe^e"* -- 4 . Then, substituting this expression in Eq. (|24[) we obtain 

iV 

= E s « eifM + ^ e_flVt ) ■ ( 26 ) 

The time independent coefficients and 6^ are determined by the initial conditions at t = for B^ v (t) and B^it). 
From Eqs. (|2"0|) and (f2"Tj) we find that these initial conditions are 

- /rT 

5^(0) = -iEVT*^- (28) 

Using the above initial conditions in Eq. (|26[) and the orthonormality property of the matrix we obtain = 
and = ^== . Replacing these values in Eq. (f26|) we get 

^(*) = E-^ e_<nr *- ( 29 ) 



r=o ^ 2 ^ r 

Using Eq. ((29)) in Eq. (|2T|) we can get easily, 

N 



^(*) = E-M*K, (30) 



i/=0 



where 

N 
r=0 

Now, we can compute the time evolution of the expectation value corresponding to the dressed occupation 
number operator h^(t) — a'^(t)a'(t), 

n„(t)=Tr[d$(t)d,' lt (t)fr i ®fo] . (32) 

where po is the dressed density operator corresponding to the atom and pp is the dressed density operator for the 
thermal radiation field and given by Eq. (|f 6[) . To compute the trace in Eq. (|32p we choose the basis |no, n%, n^)d- 
From Eq. (f3T)|) and its hermitian conjugate we have 



JV 
N 

= El/^WI a a^ + E/^,(*)/^(*)aH- ( 33 ) 

In the basis \tiq, rii, njy)d the second term in the above equation gives no contribution for Eq. (|32[) . Then, 
replacing Eq. ({53")) in Eq. ([3"2"| we obtain easily, 

iv 

«,*(*) = I/moWI^o(0) + ^ \Uk(t)\ 2 n k (0) , (34) 
where the initial distributions for the dressed atom and field modes are given respectively by 



n (0) = Tr (a$a' po) 

OO 

= d (n\p \n) o 



(35) 
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and 



n k (0) = 



Tr k (a'la' k , 



-/3^ fe (a'Ja' fe + l/2) 



Tr fe (e-^ k ^^ +1 / 2 >) 
1 



1 



(36) 



Setting p = in Eq. (|34p. we obtain for the time dependent thermal expectation value of the occupation 
number operator, corresponding to the atom, 



«o(t) = l/oo(*)| 2 "o(0) 



N 

£ 

k=l 



\fok(t)\ 2 n k (0) 



(37) 



In early references it has been showed that |/oo(i)| 2 is the probability of the atom to remain at time t in the first 
excited level, whereas |/ofc| 2 is the probability decay of the atom from the first excited level to the ground state by 
emission of a field quanta of frequency tu k [EH HOI 121] • Then, Eq. (f3"T| suggest a clear physical interpretation in 
terms of these probabilities. Also Eq. ([51)1 can be interpreted in the same way. 

For the frequency field modes given in the paragraph after Eq. ([3j and in the continuum limit L — > oo the 
coefficients foo(t) an d fok{t) are calculated in Appendix . We obtain the following values [Eqs. ([72")) and ([73)1] 



and 
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/oo(t) = (1 - ^) e --*-^/ 2 + 2igJ(t) 
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(39) 



where k = 



J(t) 



dy w 



and 
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/„ ~ a [(y 2 + w 2 ) 2 - n 2 g 2 y 2 ]{y 2 + uj 2 ) ' 
Replacing Eqs. (|38|) and (f39|) in Eq. (|34|) we obtain in the continuum limit Aoj — * 0, TV — * oo, 



where 



P6»(t) = 



Poo(t) 



2-w 2 J 



p -*gt 
2 L 



,2 p -ngt 



2gJ(t)e- 7Tgt / 2 
p-ngt/2 



2sin(Kt) + — cos(«i) + Ag 2 J 2 (t) 



(40) 



(41) 



(42) 



(43) 



kK(uS) 



K 2 (lo) 



[2k(w 2 - w 2 ) 2 + ir 2 g 2 u(uj 2 + lu 2 )] cos[(cj - 



+7rg(u; 2 — cl>q)(cl> 2 + — 2kw) sin[(cj — n)t\ + 2gl(u, t)K(ui) [2k(w 2 — cj ) sin(/ci) 



+ 7rg(w 2 + uJq) cos(«i)] ) + igK^^ [(ui 2 - uJq) sm(ujt) + irgu; cos(ivt)] + 4Kg 2 I 2 (u, t) J- , (44) 



and 



K(u) = {uj 2 - uj 2 ) 2 + TT 2 g 2 u; 2 



(45) 

Note that in the limit t — > oo Eq. (|42[) have a well defined limit, that is, the atom reaches a final equilibrium 
state. Also, in this limit the term Poo(t), proportional to no(0) vanishes, that is, the final equilibrium distribution is 
independent of the initial atom density operator, po, it depends exclusively on the thermal field degrees of freedom. 
As showed in Refs. [HI [20l EI] Poo(t) goes to zero almost exponentially in a time of the order n/g. Taking t — > oo 
in Eq. (g^J) we get 



n (oo) = 2g / dw-r—z ^5 „ „ „. , n -r 

Jo [(w 2 -cj ) 2 +7i- 2 5 2 w 2 ](e^ - 1) 



(46) 



Now, the question is about the physical meaning of the equilibrium value given by Eq. (|46[) . To answer this 
question we compute the thermal expectation value of the number operator a G Ja' , in the case in which the atom- 
electromagnetic field system is at thermal equilibrium at some given temperature In this case the density 
operator is given by 



-OH 



TT(e- 6H ) 



where H is given by Eq. ([3]). We want to compute 



n 



Tr(e- e ») 



To compute above expression we write H , as 



N ( l\ 
H = J2[ AtA r + -\Sl r , 



(47) 



(48) 



(49) 



where A r and A\ are the normal annihilation and creation operators and given by 



A r = 



1 - n r s 

1 - /?V - 

7 =P r+l] J-Qr 



(50) 
(51) 



Now, using Eq. (JTTJ> and from Eqs. (fI7 |) -([T8 | and ([50 ]) - ([5l] ) we find that 

N N 

^=E*^» ~4 = H tr A- ( 52 ) 

r=0 r=0 

Using above expressions in Eq. (|48p and computing the trace by using the basis |no,Tii, ...,njv) c , that are eigen- 
vectors of H, we find easily, 



= E (53) 

r=0 



and in the continuum limit we get [see Appendix, Eq. (|75[) 



n = 2g dx-rj-2 rj- — 2 2 2U 9:c tt ■ (54) 

In the case in which 6 = (3, Eqs. ([54"]) and (]4l)]) are identical. Then, we conclude from above calculations that the 
atom reaches a final thermal equilibrium distribution, it thermalizes with the thermal radiation field at temperature 
Note that for weak coupling g <§C ljq, we can obtain from Eq. (|46|) or ([54]), 



n (°°) ra ' ( 55 ) 



a Bose-Einstein distribution, an expected textbook result. 



4 conclusions 

In this work we have showed that an atom (approximated by the dressed harmonic oscillator) initially in any 
arbitrary state and suddenly coupled to a thermal radiation field, evolves in time to a final thermal equilibrium 
state. The mean time, necessary for this to occur, can be roughly estimated from Eqs. (j4"5]) -(|44 p and is of the order 
7r/<7, an intuitively expected result. Also, we have found a physically suggestive result for the time evolution of the 
thermal expectation value of the dressed occupation number operators, Eqs. ([M]) and Eq. ([57]) . In general, this 
time evolution is given in terms of the time dependent probabilities associated with the emission and absorption 
of field quanta. 
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A The continuum limit 

We want to compute, in the continuum limit, sums of the type 

JV 



r=0 

where TZ^ u {£i) is an analytic function of Q. For this end we define a function W(z), 

N 2 2 

W{z) = z >-u>l + Y J 4^ ■ (57) 



col — z 
k k 



From Eqs. ([7J and ((5J) we can note that the f2 r 's are the roots of w(z). For complex values of z and using 
i] 2 = 2gAui, we can write Eq. ([57]) in the continuum limit as, 

W{z) = z 2 -ul + 2gz 2 f°^% • (58) 



For complex values of z the above integral is well defined and can be evaluated easily using Cauchy theorem, 
obtaining 

W{z) = { **+W*-«§,M*)>0 (59) 
w 1 z 2 - igirz - uj z , Im(z) < . v ' 



We start by computing Roo(t), 
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Roo — y^(tp) 2 7?-oo(^ r ) 
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From the expression for £q, given in ([7]) and Eq. (|57[) it is easy to show that, 

2 _ 2fl r 

where the prime means derivative with respect to the argument. Since the f2 r 's are the roots of W(z), we can write 
Eq. (JBDJ) as 
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ic W{z) ' 

where C is a counterclockwise contour in the z-plane that encircles the real positive roots fi r , that is, a contour 
that encircles the real positive axis. The integral in Eq. (152^) can be evaluated choosing a contour that lies just 
below and above of the real positive axis. Below the real positive axis we have z — a — ie and above z — a + ie, 
where a is real positive and e — > + . Then, we have for Eq. (|62|) . 
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W(a-ie) W(a + ie) 

From Eq. (fS"9"|) we get for W(a — ie) and W(a — ie) respectively in the limit e — » + , 
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W(a + ie) 
W(a - ie) 



a 2 — ujq + igna , 
a 2 — u)q — igna . 



Taking the limit e — > + in Eq. (|6"3"|) and using Eq. we get 

a 2 TZ 00 (a) 
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(a 2 — luq) 2 + g 2 n 2 a 2 
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As a check that Eq. (|65[) is correct we take the case TZqq = 1 and using Cauchy theorem it is easy to show that 
above integral is 1, as expected from the orthonormality property of the matrix 
Next we compute J2ofe(X)> 
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Using the expressions for tg and t k , as given by ([Jj), in Eq. ((66)) we obtain 
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where in the second line the pole at z = oj k gives a zero contribution since W(u>k) as given by Eq. ([57)1 or ()58|) is 
infinity. Evaluating Eq. ([6"Tj) by choosing the same contour as in the evaluation of Rqq (t) we get 



Rok = —-r^- I da 
Using Eq. (JB^I in Eq. (|BSJ) we obtain, 
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where K = y u> 2 — \g 2 - To check the validity of Eq. ([M]) we take IZok = 1 and using Cauchy theorem it can be 
proved that the integral vanishes as expected from the orthonormality of the matrix {i^}. 
Now, it is straightforward to compute the coefficients / M „ (t) 



N 



r=0 



in the continuum limit. Taking /i = v = in Eq. ([70)) and using Eq. (|65p we get 



from which we find 



Mt)=2gl dx {x2 _ , (71) 
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/ooW = (1 - -£ )e -*«-«*l* + 2l g / y m , (« 2 > 0) . (72) 

Taking it = 0, ^ = fc in Eq. (JTUJ) and using Eq. we get 



/ofc 



<~/x 



(x - ^ _ K )(a; _ isa + K )(a, _ ie _ Wfc )(a; _ ie + Wfc ) 



(x+ ^ - k)(x + ^ + k)(x + ie- w fe )(x + ?:e + w fc ) 



(73) 



We can integrate Eq. (|73|) in the complex plane by using cauchy theorem. We choice as the closed contour of 
integration, the path that goes in the real axis from to oo, then go to the negative imaginary axis along the part 
of the circle with radius R — > oo and argument — tt/2 < 6 < and closes the contour along the imaginary axis from 
—ioo to the origin. Note that inside the contour of integration only the second term in the bracket of Eq. (f6"9")l has 
two poles at —igir/2 + n and — ie + lu^. Then, we get for Eq. (|73[) 



fok(t) = 



To compute 



— < (-2Z7T) 

ITT 



(HL£ _ ]^g-tKt-?r»t/2 







2[( K - H2)2 _ w 2] 2[wg + i7rsw fc -wg] 



(-y 2 + Trgy - W 2 )(?V 2 + W 2 ) (_ y 2 _ _ w 2)( y 2 + w 2j 



(1 — ing^ c -iKt-nQt/2 



-iiAj k t 



(/c-i^)2] - c^ 2 + i7r 5 u; fe ] 



+ 2ig?7Wfe / 



y [(y2 +CJ 2 )2 _ 7r252j/2](j/2+w 2 r 



(K 2 > 0) 



n 



N 

E 

r=0 



(t5) 2 



in the continuum limit we use Eq. (|65p . obtaining 

/•OO 

no = 2g / g?x 



[(x 2 - Wo )2 +Tr 2 g 2 x 2 }{e 9x - 1) 



(74) 



(75) 



(76) 
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